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Abstract. Motivated by brane physics, we consider the non-linear Dirac-Born-Infeld 
(DBI) extension of the Abelian-Higgs model and study the corresponding cosmic 
string configurations. The model is defined by a potential term, assumed to be of 
the mexican hat form, and a DBI action for the kinetic terms. We show that it 
is a continuous deformation of the Abelian-Higgs model, with a single deformation 
parameter depending on a dimensionless combination of the scalar coupling constant, 
the vacuum expectation value of the scalar field at infinity, and the brane tension. 
By means of numerical calculations, we investigate the profiles of the corresponding 
DBI-cosmic strings and prove that they have a core which is narrower than that of 
Abelian-Higgs strings. We also show that the corresponding action is smaller than 
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in the standard case suggesting that their formation could be favoured in brane 
models. Moreover we show that the DBI-cosmic string solutions are non-pathological 
everywhere in parameter space. Finally, in the limit in which the DBI model reduces 
to the Bogomolnyi-Prasad-Sommerfield (BPS) Abelian-Higgs model, we find that DBI 
cosmic strings are no longer BPS: rather they have positive binding energy. We thus 
argue that, when they meet, two DBI strings will not bind with the corresponding 
formation of a junction, and hence that a network of DBI strings is likely to behave as 
a network of standard cosmic strings. 

PACS numbers: 98.80.Cq, 98.70.Vc 
1. Introduction 

The (Wilkinson Microwave Anisotropy Probe) WMAP5 results [D [21 [3l H El E] give 
strong indication in favour of cosmic inflation over other mechanisms for the production 
of primordial fluctuations [7j. Since inflation generally takes place at high energy, 
recently there has been a flurry of activity in constructing models inspired by or derived 
from string theory (for recent reviews see e.g. Refs. [HI El [101 [H])- In a large category 
of these models, particularly brane-antibrane inflation [121 [131 [HI [HI [161 [13 [HI [HI [20] 
and D3/D7 inflation [211 [22], the end result of the inflationary phase is the creation of 
D-strings (as well as potentially F-strings [22]), interpreted from the four- dimensional 
point of view as cosmic strings. Since cosmic strings are strongly ruled out as the 
main originator of primordial fluctuations, the D-(and indeed F-) string tension is 
severely constrained and (under certain assumptions) such that G^n < 10~^ in order 
to preserve the features of the WMAP5 results [21]. Nevertheless, the existence and 
possible detection of the effects of D-strings in the aftermath of an era of brane inflation 
could be a testable prediction of string theory. 

D-strings themselves have been conjectured to be in correspondence with the D- 
term strings of supergravity [25]. One of their remarkable properties is that they satisfy 
a Bogomolnyi-Prasad-Sommerfield (BPS) condition, i.e. they have no binding energy 
and preserve 1/2 of the original supersymmetries. They also carry fermionic zero modes 
and are therefore vorton candidates, leading to possible interesting phenomenological 
consequences [26] . 

The identification between D-term strings and D-strings has been made in the 
low energy limit, when field gradients are small. Inspired by the case of open 
string modes which can be effectively described by a non-linear action of the Dirac- 
Born-Infeld (DBI) type, in this paper we construct models of cosmic strings which 
depart from the low energy approximation and generalise the Abelian-Higgs model 
to a non-linear one. We will call the resulting topological objects 'DBI-cosmic 
strings', and they are exact solutions of the generalised non-linear DBI action. The 
action we consider is very different from others which have been discussed in the 
literature, Refs. [23, [23 [221 [33 [3ll [32], and in particular does not lead to pathological 
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configurations. In the limit of small field gradients our DBI strings reduce to Abelian- 
Higgs strings. We construct the DBI string solutions numerically in a broad range of 
parameter space, using two numerical methods: a relaxation method and a shooting 
algorithm. In this way, we show, in particular, that DBI strings with a potential term 
corresponding to the BPS limit of the Abelian-Higgs model are no-longer BPS. More 
specifically, /i2n > 2/i„, where /i„ is the action per unit time and length for a string with a 
winding number n: the equality only holds in the low-energy limit. Borrowing language 
from the standard cosmic string literature [331 EH [35] , the strings are therefore in the 
type II regime (though the deviations from BPS are small, in a sense we will quantify). 
The network of strings produced will therefore not contain junctions, and all the strings 
will have the same tension /i„=i. In the cosmological context we therefore expect the 
DBI-string network to evolve in the standard way e.g. j36j, containing infinite strings 
and loops, radiating energy through gravitational waves, and eventually reaching a 
scaling solution. 

The paper is organised as follows. In section 12.11 we recall the properties of 
Abelian-Higgs cosmic strings, while their realisation in the D3/D7 system is discussed in 
subsection 12. 2[ In section [3l we first briefly review the different non-linear actions which 
have been put forward so far in the literature. Then in subsection 13.21 we motivate 
and present our proposed non-linear action for cosmic strings, which we expect to be 
applicable when field gradients are large. In section HI we study the DBI-cosmic strings 
solutions analytically and numerically. In subsection 14. ![ we present simple analytical 
estimates which allow us to roughly guess the form of the DBI string profiles and, 
in subsection 14. 2[ we compute them numerically by means of two different methods 
(shooting and over relaxation). In section [5] we briefly summarise our main findings and 
discuss our conclusions. Finally, the appendix gives the full non-linear structure of the 
DBI cosmic string action. 



2. Abelian-Higgs Cosmic Strings 

2.1. The Abelian-Higgs model 

We begin by recalling briefiy the properties of standard Abelian-Higgs cosmic strings, 
and at the same time introduce our notation following Ref. [3^ though we use the 
signature ( — h ++). 

The Abelian-Higgs model is governed by the action 



S 



d X^ 



-9 



(D^S)(D^S)t + -F^^F^-" + Vm) 



where the potential is given by 

A 



(1) 



(2) 



In Eq. ([T]), Dn denotes the covariant derivative defined by D^^ = — iqA^ with 
the vector potential, F^y = d^j_A„ — di,A^, and q is the gauge coupling. The potential is 
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Figure 1. Left panel: the tension of n = 1 and n — 2 Abelian-Higgs strings as 
a function of /3^^ = 2q^ /\. Right panch the binding energy /i2„ — 2/i„ for n = 1 
Abehan-Higgs strings as a function of (3^^ . When the BPS condition /3 = 1 is satisfied, 
/i„ = 2Trnr]^ so that fj,2n — 2/in = as can be verified in the figure. 



characterised by two free parameters: the couphng A > and an energy scale rj. It is 
useful to introduce the dimensionless coupling 

/J-^ = 4. (3) 

where the Higgs mass is rris = \/Xri/y/2 and the vector mass nig = qrj. 

Due to the non-trivial topology of the vacuum manifold, after gauge symmetry 
breaking the model possesses vortex (or cosmic string) solutions for which the scalar 
field can be expressed as 

S(r,^)=r7X(p)e^"^ (4) 

where we have used the cylindrical coordinates, and the cosmic string is aligned along 
the z-axis. Here n is the winding number proportional to the quantised magnetic flux 
on the string, and we have defined a rescaled radial coordinate 

p = X^/^r]r , (5) 

with X (p) —>■ 1 at infinity, while X{0) = 0. In the radial gauge, the only non- vanishing 
component of the vector potential is Ag{p) with ^^(O) = 0. We define 

Q = n- qAg , (6) 

so that the tension, defined to be the action per unit time and length p = —S/dtdz, 
can be expressed as 



Pn{(3) = ^nrj / dpp 



\dp J p^ \ dp 

= 27:r,^g^ {(3-') . (7) 
The function Qn is plotted in Fig. [1] (left panel). 
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In the BPS case, /3 = 1, the tension given in Eq. ([7]) can be re-expressed as 




(8) 




(9) 



On inserting back into Eq. ([H]) one finds 



(10) 



so that gn (1) = n. 

The functions gn=i and gn=2 are plotted in the left hand panel of Fig. [1] 

whereas the binding energy fi2 — 2yUi is shown in the right hand panel. In the BPS limit, 
gn (1) = n, and the force between vortices vanishes [37]. For P < 1, the strings are type 
I with a negative binding energy, while for f3 > 1 they are type II with positive binding 
energy. Type I string therefore attract and can form bound states or 'zippers' [38] 
linked by junctions. Zippers may form (in a certain regime of parameter space) when 
two strings in a network collide, Refs. [391 SO]- A network of type II strings, on the 
other hand, contains no junctions and the strings all have the same tension fin=i- 

We now outline how BPS Abelian-Higgs cosmic strings form in certain string theory 
models of infiation. This will motivate our discussion of DBI strings in section [3l There 
we will see that DBI strings with (3 = 1 have a positive binding energy and hence repel 
each other. 

2.2. BPS Abelian-Higgs strings in the D3-D7 system 

In string theory models, cosmic strings can form after the end of inflation [Sj [H [101 El • 
the D3/D7 system [211 [22] in particular, the two branes attract during the inflationary 
period and then eventually coalesce forming D-strings. The whole picture (inflation 
and string formation) can be described in terms of the fleld theoretical D-term hybrid 
inflation [25]. In this language the D-strings have been conjectured to be analogous to 
D-term strings, and furthermore — as we now outline — the strings are BPS Abelian- 
Higgs strings. 

In the D3/D7 system there are three complex flelds [211 122]: the inflaton (p and 
the waterfall flelds 0^. In string theory, is the interbrane distance and 0^ are 
in correspondence with the open strings between the branes. In the supersymmetric 
language, the Kahler potential is 




(11) 
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leading to canonically normalised fields. Notice that the inflation direction is invariant 
under the real shift symmetry ^ + c, thus guaranteeing the flatness of the inflaton 
direction at the classical level [H]. The superpotential is 

W = X(j)(j)+(j)- . (12) 

During inflation, the U{1) symmetry under which the waterfall flelds have charges ±1 
is not broken, i.e 0^ = 0. The scalar potential is flat and picks up a slope at the 
one loop level. This is enough to drive inflation. As (p decreases, it goes through 
a threshold after which the waterfall fleld 0"'" condenses and the inflaton vanishes. 
This corresponds to the coalescence of the D3- and D7-branes. The effective potential 
describing the condensation is given by the D-teim potential (the F-terms all vanish 
when = 0, 0~ = 0) 

Vn = '-^m'-0"- (13) 

The term ^ is called a Fayet-Iliopoulos (FI) term [42j. As 0+ condenses and (0+) = y/^, 
cosmic strings form interpolating between a vanishing fleld in the core and at 
infinity. These cosmic strings are BPS objects preserving one half of the original 
supersymmetries. Their tension is known to be /i„ = 27m^ [25]. As a consequence, 
there is no binding energy as fi2n = 2/i„. 

In fact [25], the D-term string model of D-string formation is nothing but an 
Abelian-Higgs model with particular couplings 

C = D^0+ ^ + ^F.^F^-' + Vd , (14) 

where D^cj)^ = {d^ — iAfj^cj)'^. Upon rescaling — > Af^/g and comparing with ([T]), this 
leads to the identification 

q = g, X = 2g\ i = r,\ (15) 

corresponding to /3 = 1 and, hence, a BPS system. This explains why one recovers 

Moreover, the energy scale can be given a stringy interpretation. Indeed, one 
can show that the Fayet-Iliopoulos term is related to internal fiuxes on D7-branes |43j . 
For this purpose, let us consider a ten-dimensional metric in the form 

ds^ = gf.^dX^'dX'' + gpgdXPdX'^ + gijdX'dX^ , (16) 

with /i = 0, ■ ■ ■ , 3, p = 4, ■ ■ ■ , 7 and i = 8,9. The quantity g^i, is the four- dimensional 
metric and gpg and gij the compactification six- dimensional metric. This corresponds to 
the metric on x compactifications for instance. The internal dimensions of the 
D7-brane are the coordinates a = (/x, p) while the D3-brane lies along the /j, coordinates. 
We denote by Ty the brane tension and gs the string coupling. The four-dimensional 
gauge coupling is given by 

1^^. (17) 

9^ 9s 
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where the string length is 4 = \/y2Tm' and V4 = / d'^x-^det Qpq is the volume of 
the four compact internal dimensions of the D7-brane. Consider now a dimensionless 
magnetic flux Fpq along the internal dimensions g) = 4, ■ ■ ■ , 7, of a D7-brane. Then 

^'=2^ I '^"^V^^FpqF^" ■ (18) 
Notice that the absolute value of the FI term is not fixed, it can be decomposed as 
= C/i^ where the prefactor depends on Fpq. 
In the following, we will consider cosmic string models for which the canonical 
kinetic terms have been replaced by a non-linear term of the DBI type. In effective 
actions describing string theory phenomena, and particularly brane dynamics, such a 
replacement is mandatory as soon as the gradient terms in the effective action become 
large. Indeed, the DBI action usually describes the dynamics of the open strings in 
correspondence with the brane motion (such as the 3-3 and 7-7 open strings in the 
D3/D7 system). As we have recalled, the formation of cosmic strings in the D3/D7 
system is governed by the 3-7 strings of no obvious geometric significance. In such a 
situation, and assuming that there could be higher order terms correcting the lowest 
order Lagrangian, the effect of the higher order corrections to the kinetic terms (terms 
in \D(f)'^\'^P, p > 1) would be to induce modifications of the cosmic string profile and of 
the tension. 

In the following, we do not restrict our attention to a particular setting such as the 
D3/D7 system. We discuss possible non-linear extensions of the Abelian-Higgs system 
and then motivate a specific, well defined form. We then analyse the departure from 
the BPS case induced by the higher order terms. 

3. DBI Cosmic Strings 

3.1. Non-standard actions for cosmic strings 

As discussed in the previous section, we are interested in situations in which gauged 
cosmic strings form when higher order corrections to the kinetic terms in the action 
cannot be neglected. In the absence of an explicit derivation from, say, string theory, 
we take a phenomenological approach (which, however, is strongly inspired by string 
theory). This is presented in subsection 13.21 We will construct an action [given in 
Eq. ( l35|) and reproduced below] which satisfies the following two criteria: Firstly, the 
Abelian-Higgs limit should be recovered when gradients are small. In particular, the 
action should be a continuous deformation of the Abehan-Higgs model. Secondly, the 
resulting cosmic string solutions should have no pathological and/or singular behaviour 
as the model becomes more and more non-linear (in field gradients). In the remainder 
of this subsection, we compare our action, Eq. ( 14T1) : 

Sex I d^xj^-det [g^, + (D(^0)(D,)0)t + F^,] - ^ 
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to other non-linear actions which have been considered in the hterature and which do 
not satisfy the above criteria. 

In Refs. [271 EHl SO] an attempt to construct a non-hnear model for (electrically) 
charged vortices in (2+1) dimensions uses an hybrid approach with a (truncated) Born- 
Infeld action for the gauge field, a standard linear action for the Higgs field, and a 
Chern-Simons tern|j|: 



S 



1 



cr 



1 -I F Pm^^ — 1 



K 



e^'''PA F 



+ l\Dct>\' + V{4>) 



(20) 



where a is a parameter of dimension two and |-D0p = (-0^0) (-0^0)^ At a threshold 
a = o"c corresponding to the very non-linear regime, the gauge field becomes singular at 
the origin of the vortex whilst, below the threshold, no solution exists. Incorporating 
the Higgs kinetic terms into the square root, while dropping the Chern-Simons term, 
leads to the following expression 



S= j dSv^ja^ + i^F^uF^"" + ^ \D<P\^ + ^^F,F,D^<P'^D\ 



2a^ 



1/2 



V(0) 



(21) 



where F^ = e^^ypF^f j1. With such a very particular form (which differs from the one we 
will propose shortly), one does not find finite energy solutions. 

A model for (global) cosmic strings was proposed in Ref. [29] with an action given 

by 



S 



d^X\J—g 



1 + V{^) 



(22) 



(this is identical to Eq. (I2T1) in the global limit). For this model the solutions become 
multi-valued and undefined at the origin as soon as the magnitude of V (0) becomes 
sufficiently large [29] . 

In view of these negative approaches where singularities and pathologies abound, 
one could be tempted to think that non-linear cosmic string actions all lead to these 
problems. Fortunately, a well-behaved action has been suggested by Sen [H] and studied 
in Ref. [15| in the case of D-strings obtained at the end of D-D inflation. In such a system, 
hybrid inflation occurs and the role of the waterfall field is played by the open string 
tachyon T with a charge ±1 under the UiX) gauge groups of the D3- (respectively D3-) 
brane. When the two branes coincide the effective action reads 



S = -tJ d^a:^ (T,Tt) [^det {-g+) + ^det {-g 



(23) 



I No charged vortex solutions exist when the Chern-Simons term is absent. 
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where g^^ = g^^ ± (-iF^y + {^D^TDyT^ + D^T'^DyT^ /2 and V is the tachyon potential. 
Since det(— (7+) = det{—g~) (see the Appendix), the action reduces to 



S = -2T3 j d^xV (T, T^) 



-9^) , (24) 



and it admits BPS vortex solutions. 

Topological defects in models with general non-linear kinetic terms were studied in 
Refs. [3T1[32]. The action proposed in pT] for global topological defects differs from the 
tachyon action above (in the global limit) as the potential is added to the generalised 
kinetic terms: 



S= / d^x 

Here 



(25) 



X = -{d^cPad^c^a) (26) 



1 

2' 

is the standard kinetic term, K{X) is some non-linear function, 0^ is a set of scalar fields, 
M has dimensions of mass, and the potential term provides the symmetry breaking term. 
One of the main restrictions imposed in |3T] on the form of the non-linear function K{X) 
is that K{X) should have a canonical asymptotic form, K{X) ~ X as X ^ 0. However, 
for large gradients K{X) could deviate considerably from the canonical kinetic terms. 
The former requirement implies a non-pathological behaviour of solutions far from the 
defect core, while the different possibilities for K{X) at infinity leads to deviations of the 
defect from the standard case inside the core. The action fl25|) leads to non-pathological 
solutions for so-called fc-defects — domain walls, cosmic strings and monopoles — whose 
properties can differ considerably from those of standard defects. 

A gauge version of action ( |25l) was considered in Ref. [32]: for a complex scalar 

field 



S= d^x 



(27) 



with 



It was shown that non-pathological cosmic string solutions exist at least for some choices 
of the non-linear function K{X) [32] . 

In the following we will motivate and study a non-linear extension of Abelian-Higgs 
model which retains some of the properties of the tachyon and fc-defect models. The 
potential will be additive as in the fc-defect case while the kinetic terms have a DBI 
form as in the tachyon case. However, the kinetic terms are not a function solely of 
X anymore: they differentiate between the radial and angular gradients of the defects. 
This springs from the origin of the kinetic terms as induced from the normal motion of 
a D3 brane embedded in a larger space-time. 
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3.2. A DBI Action for Cosmic Strings 

We now turn to the action that we propose in this article. Consider a brane model in 
which cosmic strings appear as deformations of a brane. (In a sense the brane becomes 
curved with a puncture at the location of the string, as we discuss.) To do so, consider a 
ten-dimensional setting as is natural for brane models derived from or inspired by string 
theory. We choose a non-warped compactification and write the ten-dimensional metric 
in cylindrical form 

ds?o = g^dX^dX'^ = dsl + 2g^-pdZ''dZ~^ , (29) 

where 

dsl = gt^udX^dX"" = -(dX°)2 + dR^ + RMQ^ + dZ^ . (30) 

The metric along the internal dimensions g^f^ {a = 5, 6, 7) is kept arbitrary, 
i.e. Hermitian and positive definite, and we have assumed that the six-dimensional 
manifold is complex (it could be a Calabi-Yau manifold) therefore having complex 
coordinates. The complex coordinates are crucial to analyse cosmic strings. 

Consider the DBI action for a three-brane embedded along the first four coordinates 



S = -T j d^x^-det {~g,u + iir,,) - j d'^x^V [VtZ' 



(31) 



where T is the brane tension, jF^j, is the field strength on the brane (and has dimension 
two), distances have dimension minus one and has dimension one. We have included 
a potential for the deformations of the normal directions to the three-branes. As 
suitable when the normal directions are charged under the world-volume gauge group [in 
this case the local U{1) on the brane], we include a covariant derivative in the definition 
of the induced metric 

~g,. = + 9a(3 (V^Z'^V^Z'^ + V^Z^V^Z^') (32) 

with 

= d^- iqAf, . (33) 

Clearly, when the gauge fields vanish, g^^i, is simply the induced metric on the brane. A 
similar extension of the induced metric to charged fields has already been introduced in 
the context of N-coinciding D-branes [16] with the corresponding non-Abelian SU (N) 
gauge theory. There the brane coordinates are in the adjoint representation and have 
kinetic terms involving the SU (N) covariant derivative [16] . We extend this procedure 
to the DBI cosmic string situation with a f/(l) gauge group|§| 

When the six-dimensional metric is nearly flat g^f^ = locally, the action becomes 

S= - T j d'a:| y/- det [g^, + {V^Z-V.Z^ + V^Z^V.Z^) + £2^^^] 
- v^} - j d'x^v{VTZ") , 



(34) 

§ In the D-brane context, the brane fields do not carry any U{1) charge as they belong to the adjoint 
representation. 
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where, as usual, we have subtracted the action of the "flat" brane so that the Abehan- 
Higgs model is recovered when gradients are small. In the following we suppose that 
only one normal direction is excited and deflne ^ = . The resulting action is given 

by 

S= - T j d^xj ^ - det [g,, + (I)(^<l>)(I?,)$)t + £2^^^] _ ^ 



T 

Notice that, in the above equation, the potential [i.e. V{x) as a function of x\ is still 
given by the expression ([2]) and, therefore, contains the parameter A. When the complex 
scalar fleld $ vanishes, Eq. ( l35l) describes Born- Infeld electrodynamics |47j . 

We now discuss action (135|) in detail. In particular we compare its properties to 
those of the Abelian-Higgs action ([T]) discussed in section O and then we construct the 
static cosmic string solutions of the action. 

A flrst important property of Eq. ( l35i) is that, to leading order in derivatives, it 
reduces to the standard action ([1]) on identifying 

S = Vt* , (36) 

and redeflning the charge according to the following expression 

9 = (37) 

together with the gauge fleld 

Hence, if the spatial derivatives characterising DBI-strings are small (we shall discuss 
whether or not this is the case below), their properties should to be very similar to 
Abelian Higgs strings. More generally, however, and as discussed in detail in the 
Appendix where we calculate the determinant explicitly, Eq. ( l35l) contains terms of 
higher order in covariant derivatives as well as numerous different mixing terms between 
T and V (suitably contracted). These extra terms could signiflcantly change the string 
solution and the resulting strings' properties relative to the Abelian Higgs case. It 
follows from this that our action is very different from that considered by Sarangi in 
Ref. [29], even in the global case. As a consequence we will flnd non-pathological cosmic 
strings solutions with a continuous limit to Abelian-Higgs strings. 

We now focus on the cosmic string solutions of Eq. (!35|) . For this purpose, flrst it 
is useful to pass to dimensionless variables, denoted with a hat. Explicitly we set 

^ = ^.u^llT,., V.^IJ),, ^^7f7' ^^^^ 



= ^ . (3^ 



as well as 



vm^^l^'-v'Y , (40) 
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where A = AT£g, so that the action becomes 



S 



Tit 



d'^xi ^ -det 



+ v^\/(|<l|)| 



(41) 



We now follow the procedure outlined in section 12.11 for Abelian-Higgs cosmic strings, 
however for action fHTj) . In dimensionless cylindrical coordinates, ds^ = —dP + df^ + 
f^d^^ + dz"^ and in the radial gauge {Af = 0), the cosmic string profile is 

$ = fiX{p)e'^' , Q{p) =n-qAe (r) , (42) 
where we have defined a new radial coordinate p by the following expression 

p = X^/^fir (43) 
which should be compared to Eq. ([5]). The boundary conditions on the fields are 



lim X = , lim Q = n , lim X = 1 , lim Q = . 

p— »0 p— >0 p— >oo p— >oo 

Substituting Eqs. ( H2l) and ( H3i) into ( HTl) as well as using ( l39l) . we find tha 



(-det5f^,,)7 ^ = - det 
where we have defined the 7 factor by 



fdxV 



aQ2x2\ a ^dQ 



p2 V "^P 



(44) 



(45) 



(46) 



Hence the string tension defined as —S/i1dzdt is given by 



IJ'n 



a 



dpp 



\ 



p2 p2 



a 



1 + -(X^ - 1)^ 



(47) 



where 

a = 2A?7^ , 
and, as in the Abelian-Higgs case, 

^ 2g2 2g2 



(4^ 



(49) 



Eq. (1471) is the main result of this section and represents the non-linear DBI 
generalisation of the linear Abelian-Higgs model: it should be compared to Eq. ([7]). 
Notice that it involves the single additional parameter a which measures the deformation 
from the Abelian-Higgs model, since Eq. f HTl) reduces to the tension of Abelian Higgs 



Note that 7^=1? where T) is discussed in the Appendix. 



Dirac Born Infeld (DBI) Cosmic Strings 



13 



strings in the linear limit a 0. As discussed in section I2.H Abelian-Higgs strings 
are BPS when (3 = 1, and hence for (3=1, DBI-cosmic strings with tension given by 
Eq. ( H7I) are a continuous deformation of the BPS Abelian-Higgs strings. This property 
is, of course, very important and constitutes an additional motivation for the action 
given in Eq. (HTl) . 

Finally, we note that the argument of the cosmic string profile p is also identical 
to its counterpart in the Abelian-Higgs case, whatever the value of a. Hence we will be 
able to find continuous deformations of the cosmic string profiles parameterised by a 
and depending on the universal variable p. 



4. DBI String Solutions 



4-1. Analytical Estimates 

Having established the model and its action, we now turn to the solutions of the 
equations of motion. The DBI cosmic string equations follow from Eq. (I47l) and read 

2" 



Jp 



P7 1 + 



7dg 

p dp 



(3p 



dX 



P 



(X2-l)X + 7 



P 



l + a 



dxy 



(50) 
(51) 



for the scalar field and gauge fields, respectively, where 7 is defined in Eq. (H6l) . In the 
Abelian Higgs limit, a = 0, Eqs. (|3U]) and ([SID reduce to the standard cosmic string 
field equations for which 7 = 1. Deviations from Abelian Higgs strings will occur if 
7 < 1. Notice that here the fields are purely space-dependent. For time-dependent fields, 
and particularly in DBI inflationary cosmology with infiaton 0(t) whose dynamics is 
described by action (HTj) in the global limit, then 7 is a generalisation of the cosmo logical 
Lorentz factor. Indeed, as can be seen from Eq. (H6l) in the case when g^^ describes 
an homogeneous and isotropic manifold, 7"^ = 1 — (j)'^/T{(j)) where T[((>) is related to 
the metric of the extra-dimensions. The difference in sign between spatial and temporal 
derivatives is responsible for the fact that deviations from standard cosmology (7 = 1) 
occur here when 7 —>■ +00 rather than 7^1. 

Unfortunately, as is clear from Eqs. (1501) and (ISTl) . the DBI cosmic string equations 
cannot be solved exactly. We have therefore carried out a full numerical integration 
of the equations of motion. For convenience, we will focus on the deformed BPS case 
where (3 = 1 and a ^ (though (3^1 and a ^ can also been done with the numerical 
methods used here). 

Before discussing the numerical results, we analyse the asymptotic behaviour of the 
fields in order to obtain a rough understanding of the solution. We will consider the two 
limits p ^ and p — 00 and will address two issues. The first one is the non-existence 
of singularities in the core of the cosmic string. The second one will be the determination 
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of the shape of the cosmic string profile both at the origin and at infinity. In particular 
we will find that the functional form of the string profile is similar to the Abelian-Higgs 
case inside the core, the only difference springing from a-dependent factors. 

Consider first the p — >^ limit and let us examine the possibility of singular DBI 
strings deep in the string core. As already discussed, the DBI features of the solutions 
depend on 7. In particular, extreme deviations from the Abelian-Higgs case would 
appear if 7 ^ at the origin. This can only happen if the derivative of X and/or Q 
become extremely large, i.e. the string becomes singular. Let us first assume that the 
gradient of Q becomes large and dominates the 7 factor, i.e. 7 ~ pdpj {^fadQ). The 
gauge equation becomes non-sensical as the left-hand of (ISTll vanishes and the right- 
hand side does not. Hence there is no regime where the gradient Q is arbitrarily large 
leading to 7 —> at the origin. We now examine the possibility that X becomes singular 
at the origin with a large gradient. In this limit, we find 

'' (52) 



aX (dX/dp) 



^ i i i p- 



2a(dX/dp)' 2n2aX2_ 

where Q ^ n close to the origin and we have expanded 7 in 1/ [a [dX/dpf] < 1 and 
p^/ (aX^) <^ 1, this last condition being the only one compatible with the condition 
on the derivative of X. Working to first order in these two parameters, the profile is 
determined by 



d 



X 



P' 



dX 



_a (dX/dp)' aX2 



' (53, 



n^aX^ a{dX/dpf_ 

Notice that to zeroth order in the two small parameters, the equation is tautological. In 
the limit p ^ with an ansatz X ~ p'' the equation of motion is satisfied for 5^ = n^. 
The only solution satisfying X(0) = is obtained for 6 = n which has finite derivative 
at the origin. This contradicts our premises and, as a result, we conclude that singular 
DBI strings do not exist. 

Having shown that the strings are not singular, we will now show that the functional 
form of the solutions is similar to the ones in the Abelian-Higgs case. Let us assume 
that, in the limit p — > 0, the DBI solutions are of the form 

X(p) = A,3,/, Qip) = n - B,,,p^ , (54) 

where p and q are two constants which we will determine below, while A^bi ^dbi 
two constants to be obtained by numerical integration; and we have taken into account 
the boundary conditions at p = 0: X(0) = and Q{0) = n. By direct substitution of 
the asymptotic form fl5^ into the equations of motion Eqs. fl50l) . flSTl) and taking the 
limit p — > 0, one can check that the correct asymptotic form for X and Q reads, 

X(p) = A,3,p", Qip)=n~B,,,p\ (55) 

Thus p = n and q = 2 and, as guessed above, the only difference between DBI cosmic 
strings and Abelian-Higgs cosmic strings close to the origin is in the numerical values 
of the prefactors A^^j and -B^^j which are a-dependent. In particular, these coefficients 
become large for large a implying that away from the origin but for reasonable and 
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finite values the gamma factor becomes noticeably different from one, i.e. the cosmic 
strings are in a mild DBI regime. 

From Eqs. f H6l) and ( l55i) it immediately follows that 7 is always finite at p = 0. 
This a salient point as it confirms that the cosmic strings constructed with a DBI action 
are non-singular at the origin. This is of course yet another argument supporting the 
fact that Eq. P7|) represents the natural DBI generalisation for cosmic strings. 

More precisely we find that as p ^ there are two possible regimes: the standard 
regime where |1 — 7I ^0 and a mild DBI regime where 7 < 1 but finite. Let us first 
analyse the global string, for which Q = 0. It is clear from Eqs. (H6|) and (155!) that 
for n > 2 the mild DBI regime cannot be realised around p 0. Note, however, that 
for a ^ 1 we find numerically that ^dbi ^ 1' which implies that away from the origin 
the gradient dX/dp becomes large, so that the solution is in the mild DBI regime. For 
n = 1 the mild DBI regime is valid starting from p = 0, if a ^ 1. For a <^ 1 the regime 
is always of non-DBI type, independently of n. 

In the case of gauge strings, the situation is similar in the limit p — > 0. Again 
for a ^ 1 the non-DBI regime is realised. For a ^ 1 we find numerically that A^^^^ 
and in Eq. ( I55l) are large. Thus the gradient dQ/dp is large too, while the terms 
proportional to dX/dp and to Q^X^ are large only for n = 1, and they are small in a 
small region around p = for n > 2. However, these terms become large away from the 
origin, since a large constant A^^^^ implies that dX/dp becomes large at some point. In 
conclusion, we find that for a large enough, the cosmic strings are in a mild DBI regime 
for finite values of p. This is confirmed numerically. 

Finally let us notice that at infinity, independently of a, both the gradients dX/dp 
and dQ/dp are small, and the cosmic string matches the standard behaviour. This is in 
agreement with the general findings for topological defects with a non-canonical kinetic 
term. 

In summary, the difference between the Abelian-Higgs and DBI strings will be small 
very far from the core of the string, while the DBI string can differ from the Abelian- 
Higgs one inside the core of the string: the larger a the larger the difference inside the 
core. 

4-2. Numerical Solutions 

As mentioned above, the equations of motion ( 150|) and ( l5Ti) cannot be solved analytically. 
For this reason, we now turn to a full numerical integration. 

As is well-known, the numerical integration is non-trivial because the boundary 
conditions are not fixed at the same point. The solutions discussed is this article have 
been obtained by means of two independent methods: a relaxation method [HI HQ], [50] 
and a shooting method. More precisely, the former is in fact the over relaxation method. 
The over relaxation method differs from the relaxation method (also known as the 
Newton iteration method) by the fact that the Newtonian iteration step is multiplied 
by a factor of 00. In the standard case, convergence for the over relaxation method is 
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Figure 2. Left panel: the profiles of global DBI strings for n — \ and for various 
values of the parameter a defined in Eq. (|48)) . Right panel: same as left panel but 
with n = 2 



guaranteed provided the over relaxation parameter uo <2 and, therefore, a good choice 
is for instance uj ~ 1.99. Here, the highly non-linear nature of the equation of motions 
may render the over relaxation method unstable. To deal with this problem, we have 
considered a "step-dependent" over relaxation parameter uj interpolating from ~ 1, 
close to the origin, and to u; — 1.3 at "infinity". As already mentioned the choice 
UJ = 1.3 -C 1.99 is due to the highly non linear behaviour of the equations. We have 
observed very severe instabilities for higher values of a;. On the other hand, the shooting 
method can be directly implemented in its standard formulation in the case of global 
strings, since there is only one integration constant to be obtained, ^dbi- While in the 
gauge case the presence of two "shooting" constants, A^^j and B^^^, makes the direct 
implementation of the standard scheme impossible, we have thus modified the shooting 
method appropriately. All in all, the two different numerical procedures, a relaxation 
and a shooting method, give the same numerical solutions, up to small numerical errors. 

Numerical integration of the equations of motion (!50|) and (!5T|) are presented and 
discussed below. 

Firstly, in Figs. [2l we consider global DBI strings (that is to say without the gauge 
field) for, respectively, winding numbers n = 1 (left panel) and n = 2 (right panel) 
and different a's. This figure confirms the qualitative statements made in the previous 
subsection. We notice that, even for "non-perturbative" values of a, i.e. a > 1, the 
difference between the standard and the DBI profiles remains quite small. Moreover, 
as announced, the maximum difference lies at a (dimensionless) radius p of order one, 
namely half way from the origin and the region where X ^ 1. Another remark is that 
the DBI profiles are always above the standard profiles. This is of course expected since 
the DBI regime means larger derivatives which, in the present context, implies the above 
mentioned property. Finally, one can check that the asymptotic behaviours discussed 
in the previous subsection are clearly observed in Figs. [2j Indeed, for n = 1, we notice 
that X{p) ~ A^BiP where A^^j is clearly a function of a (see in particular the zoom in 
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Figure 3. Left panel: the solid lines represent the profiles of the scalar and gauge 
fields of a DBI strings with n = 1 and a = 1 while the dashed lines are the profiles 
of the scalar and gauge fields of a standard Abelian-Higgs string (i.e. a — Q). Right 
panel: same as left panel but with n = 2. Notice that, on the y-axis, we have used the 
notation a = qAg/n. 

the left panel). The same remark applies for n = 2, where X{p) ~ p^. 

Secondly, in Figs. [31 we display the profiles for DBI local strings, i.e. for the scalar 
field and the gauge field, in the case where a = 1, n = 1 (left panel) and n = 2 (right 
panel). The same remarks as before apply. In particular, the scalar field profile always 
lies above its Abelian-Higgs counter part and, on the contrary, the DBI gauge field 
profile always lies inside the standard profile. As already discussed, this is because, in 
the DBI regime, the gradients are, by definition, larger than in the standard case. This 
means that a DBI string has a core smaller than an Abelian-Higgs string. As before, 
one can also check that the asymptotic behaviours are those discussed in the previous 
subsection. This is true in particular for the gauge field for which we always see that 
Q ~ n — at the origin. 

Thirdly, additional information on the profiles can be gained from Figs. HI In the 
left panel, we have compared the local and global profiles. One can notice that the global 
profile is less concentrated than the local one. Another remark is that the difference 
between the Abelian-Higgs and DBI profiles is more important in the local case than 
in the global one. In the right panel, we have compared the slopes at the origin. In 
the standard case, one has X{p) ~ ^standardP" been argued before that in the 

DBI case, one also has X{p) ~ A^^jp". We have represented the ratio ^DBiMstandard ^'^^ 
various values of a and n. One notices that the larger a, the steeper the DBI slope, 
which seems natural since the value of the parameter a controls how important the 
DBI effects are. We also remark that the same trend is observed when one increases n 
rather than a. In conclusion, from these two figures, one confirms that the deeper one 
penetrates into the DBI regime, the narrower the core of a cosmic string is. The effect 
is larger in the local case than in the global one and for large winding numbers than for 
small ones. 
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Figure 4. Left panel: comparison of the cosmic string profiles for global and local 
DBI strings with n = 1 and a = 1. The solid line represents the scalar field profile for 
a local Abelian-Higgs string while the dotted line is the corresponding DBI profile. On 
the other hand, the dashed line represents the scalar field profile for an Abelian-Higgs 
global string whereas the dotted-dashed line is the corresponding DBI profile still in 
the global case. Right panel: ratio ^DBi/^Hta„dard (^^e the text) as a function of the 
parameter a for different values of the winding number n. 



Fourthly, given the numerical solutions presented above it is straightforward to 
calculate their tension which, from the action given by Eq. fHTj) . takes the form 

/i„(X,Q) =27rr7Vn(«), (56) 

where /n(0) = n in the Abelian-Higgs case. In Fig. O (left panel), we plot the universal 
functions /„(«) for the DBI local strings. We notice that the DBI action is, for any n 
and/or a, smaller than the corresponding standard action. Moreover, at a fixed value 
of a, the (absolute) difference between the DBI and Abelian-Higgs actions increases 
with the winding number. The fact that the DBI action is smaller than the standard 
one is likely to have important physical consequences, in particular with regards to 
the formation of DBI cosmic strings. Indeed, if their energy is smaller than in the 
standard case, one can legitimely expect their formation to be favoured as compared to 
the Abelian-Higgs case. 

In the right panel in Fig. [Sj we have represented the DBI string binding energy 
/i2„ — 2/i„ as a function of the parameter a for different values of the winding number 
n. We observe that this quantity is always positive but small in comparison to one. 
Moreover, as expected since one has (/i2n — 2/i„) (a = 0) = 0, it increases with a. We 
conclude that when a 7^ 0, the DBI cosmic string is no longer a BPS object. The fact 
that /i2n > 2/i„ means that, when they meet, two DBI strings will not constitute a 
new single string with winding number 2n since this appears to be disfavoured from 
the energy point of view. This has important consequences for cosmology since the 
above discussion implies that the behaviour of a network of DBI cosmic strings will be 
similar to the behaviour of a network of Abelian-Higgs strings. This means that the 
cosmological constraints derived, for instance in Refs. [21], also apply to the present 
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Figure 5. Left panel: the solid lines represent the DBI string tension as a function 
of a for different values of the winding number n (from rt = 1 to n = 4 going from 
the bottom to the top of the plot). The dashed lines corresponds to the Abelian-Higgs 
tension, namely /i„ = 2i:Tf"n, and are thus horizontal lines located at the y-coordinate 
n. Right panel: the DBI string binding energy ^2n — 2/i„ for various n as a function 
of the parameter a. 




Figure 6. Left panel: the energy density for DBI strings as a function of p for n = 1 
for different values of the parameter a. Right panel: same as right panel but for n — 2. 



case. 

Finally, in Fig. El we represent the energy density as a function of the dimensionless 
radial coordinate p for n = 1 (left panel) and n = 2 (right panel) for different values 
of the parameter a. We notice that the DBI energy densities are usually more peaked 
than the Abelian-Higgs ones. Moreover, the larger a, the more peaked the distributions. 
The case n = 2 is particularly interesting. One observes that, as a increases, the peaks 
of the distribution are displaced towards the left, i.e. towards smaller values of p. This 
is probably due to the fact that, as discussed at the beginning of this subsection, the 
difference between the DBI and Abelian-Higgs profiles is maximum for intermediate 
values of p. 
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5. Conclusions 

We have considered a natural DBI generalisation of the Abelian-Higgs model whereby 
the kinetic terms of the Higgs fields do not lead to a hnear differential operator in 
the equations of motion. The particular form of the action is motivated by a specific 
extra-dimensional model where the Higgs field becomes a complex direction normal to a 
D3-brane. Although this model leads to nice cosmic string properties, it is not directly 
related to a string theory model. As such, the closest model of string theory which 
could have lead to such a DBI action is the D3/D7 system where BPS cosmic strings 
are formed at the end of an hybrid-like inflation phase. Unfortunately, the charged 
fields associated to the open string joining the D3- and D7-branes have no obvious 
geometric meaning and therefore do not lead to our DBI action. It would certainly be 
very interesting to see if our construction can be embedded within string theory. 

As a four-dimensional model of non-canonical type, the DBI model of cosmic 
strings does not suffer from any pathology such as divergences or non-single- valuedness 
of the field profiles (typical of other non-linear actions which have been proposed in 
the literature). Indeed we find that DBI strings can be continuously deformed to 
their Abclian-Higgs analogue. In fact, the main difference from the Abclian-Higgs case 
appears in the BPS case where the DBI strings show a small departure from the BPS 
property. In particular, we find that the string tension is reduced, a property which 
may have some phenomenological significance in order to relax the bound on the string 
tension coming from Cosmic Microwave Background (CMB) data. Moreover we find that 
the DBI strings have a positive binding energy implying that the string coalescence is 
energetically disfavoured, leading to the likely formation of networks with singly- wound 
strings and statistical properties akin to the usual Abelian-Higgs ones. 

In the present article, we have not tackled some important aspects of DBI string 
dynamics such as string scattering (for which we expect that the higher order terms 
discussed in the Appendix may play an important role), gravitational back reaction, and 
the coupling to fermions and their zero modes. This is currently under investigation. 

In summary, we have introduced DBI cosmic strings as non-singular solutions 
derived from a non-linear Lagrangian. We have studied the solutions numerically and 
found that they differ significantly from their Abelian-Higgs analogues. However, the 
network properties of these strings is almost certainly similar to those of type II Abelian- 
Higgs cosmic strings. 
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6. Appendix 

In this appendix, we study in more detail the general form of the action (l35l) considered 
in this article. Eq. fl35|) reads 

S= j d^a; I y/- det [g^, + (P(^$)(P,)$)t + P^J^^,] 

V (Vt|$|) 



T 



T j d^xv^^ [VV-l^ + 



VWT\^) 
T 



(57) 



where P is defined by 

V = det [5^, + (P'^<l>)(I?,<l>)t + + ^lJ^^'^] , (58) 

T has dimensions of (energy)"^ and 4 is a length scale. As before, = d^ — iqA^. 

Our goal is to compute and simplify Eq. fISS]) for V. As it is clear from its definition, 
this will allow us to derive a more compact formula for our action in the general case. 
In Eq. ( H71) we have evaluated the action (1571) for a cylindrically symmetric static string 
profile. In this case it takes a simple form. However, when there is time dependence 
and less symmetry — as occurs for example in string scattering — it is important to 
know the general form of the action. 

First define the following quantities 

Ar'^ = I?M$^ 5^ = A^A^ + A'^A^ , i?^^ = ^^ + J^^, (59) 

where are bar denotes complex conjugation and we set = 1 in this appendix. Note 
that by definition S^^ is a symmetric matrix and J^^'^ is antisymmetric, while S^y and 
T^v are in general neither symmetric nor antisymmetric. Denote by S the matrix with 
components S^i,^ while is the matrix with components J^^^. For integer n and p 

tr (5PJ^2"+^) = . (60) 

On the other hand, we also have 

V = det {6^"^ + R\) (61) 

x(5'^3^3+i?"3^3)(5'^V4 + ^"*/34) (62) 
which, on using the identity 

e.,a.«3«4^°^-"^^^"^-^* = - (4 - j)!j!5i^C/:::« > (es) 

gives 

I? = 1 + + R^'^o.R^'^p + R^°aR^pR''\ + R^^'^R^R^^R^^s ■ (64) 
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We now evaluate each term in the above equation. For the first (hnear in i?), it 
follows from Eqs. fl59l) and fl60l) that 

= = 2N^N'' = 2 (1"^$) (V^^)^ . (65) 

The quadratic term is given by 

= + 25l°„^^l^ + , (66) 

= i[tr^(5)-tr(5^)]-ltr(^^) , (67) 

= {N^N'^Y - (^"^") i^pN") - ^tr (^2) ^ (68) 

where to get from Eq. (1671) to Eq. (1681) we have used Eq. ( |59l) . Notice that these terms 
are compatible with the U{1) invariance of the action. The next step is to calculate the 
cubic term. It is given by 

+ J^["„J^^/3J^^^ . (69) 

The term in S"^ vanishes for the single complex scalar field studied here since, on using 
Eq. ( l59l) . it contains the contraction of an antisymmetric tensor with a symmetric one. 
Similarly S^^'^S^T^^ = = T^"" ^.Tf^ on using Eq. (160|) . Therefore, the cubic term 
takes the form 

R^''aR%R''\ = ^S^'^aO'^r^ = ^ [-ti (S) tr (^2) _^ 2tr (SJ^^)] . (70) 

Finally, the quartic term can be expressed as 

R^''o.R^pR\R^\ = S^''^Sf^(3S\S^^s + iS^'^aS'^pS^J^^^s + ^S^'' o.T%T\T^\ 
+ 6^["„5^^^^^^'^]5 + (71) 
= ^S^'^ ^S%rL,T^^ s + T'^''o.T%r<^T^^i , (72) 

since the terms on the first line in the above equations vanish, on using the same 
arguments as above. Also 

S^''^S%T\T^\ = ^|4tr (S) tr (ST^) - 4tr (j^'^S^) - 2tr (TSJ^S) 



+ [tr {S') - tr^ (S)] tr {j^') | . 



(73) 



J^^'^aJ^^J^^J^^^s = ^ [-6tr (J^^) + 3tr2 (j^^)] (74) 
Therefore, in the end, one obtains the following expression for V 
V = l + tr{S)- Itr (^2) + 1 [tr^ (^^) - 2tr (^^)] 

+ i [tr^ (5) - tr (5^)] + ^ [2tr [SJ^') - tr (S) tr [j^') 

+ i [tr (52) - tr^ (S)] tr {j^^) + tr (S) tr {SJ^^) - tr {j^^S^) (75) 
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The three terms of the first hne in Eq. fl75l) . when substituted in Eq. fl57p and on 
expanding the square-root, give the standard Abelian-Higgs model. The last two terms 
of Eq. (1751) are the standard terms of Born-Infeld electro-dynamics. Finally, as discussed 
in the main text, the factor T) and, hence, our action defined by Eq. (135|) . contains 
terms higher order in covariant derivatives as well as mixing terms between J-"^ and the 
covariant derivatives. 
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